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Significant Findings 


This paper presents one of the first algorithms to estimate snow rate over land. This 
algorithm is unique with respect to the two other early empirical snowfall algorithms in 
that it is physically-based and relates hydrometeor properties to satellite-based passive 
millimeter-wave observations. This work relies on the Advanced Microwave Sounding 
Unit (AMSU-B) millimeter-wave (89, 150, 183±1, ±3, and ±7 GHz) frequencies to 
observe the snow in the atmosphere. The higher frequencies are often advantageous 
because not only are they sensitive to frozen particles but they are also very sensitive to 
water vapor that obscures the highly variable underlying surface. In addition to snowfall 
rate, relative humidity profiles and fraction of ground covered by snow are also retrieved. 

We found that the theory of scattering by randomly oriented dry snow particles at high 
microwave frequencies appears to be better described by regarding snow as a 
concatenation of “equivalent” solid ice spheres rather than as a sphere with the effective 
dielectric constant of an air-ice mixture. An equivalent sphere snow scattering model was 
validated against high frequency attenuation measurements. Satellite-based AMSU-B 
high frequency observations during the March 5-6, 2001 New England blizzard were 
used to retrieve snowfall over land. Vertical distributions of snow, temperature and 
relative humidity profiles were derived from the Pennsylvania State University-National 
Center for Atmospheric Research (PSU-NCAR) fifth-generation Mesoscale Model 
(MM5). Those data were applied and modified in a radiative transfer model that derived 
brightness temperatures consistent with the AMSU-B observations. The retrieved 
snowfall distribution was validated with radar reflectivity measurements obtained from 
the National Oceanic and Atmospheric Administration (NOAA) National Weather 
Service (NWS) ground-based radar network. 
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Popular Summary 


Although most global precipitation occurs as rainfall, snowfall plays a significant role in 
northern latitude hydrological cycles. One of the most important challenges for the future 
is to detect snowstorms from space. This study developed one of the first non-oceanic 
snowstorm estimation algorithms. Over land, the challenges are especially great: (1) 
satellite observations used for rainfall estimation are rather insensitive to frozen particles, 
(2) land surface characteristics can contaminate the observations, and (3) the relationships 
between the physical properties of frozen hydrometeors and passive brightness 
temperatures are not well understood. This work relies on millimeter-wave frequencies to 
observe the snow in the atmosphere. These frequencies are often advantageous because 
not only are they sensitive to frozen particles but they are also very sensitive to water 
vapor (almost always present in the atmosphere). Thus, they obscure the underlying 
surface by taking advantage of water vapor screening. This study presents a physical 
model of radiation from frozen particles at millimeter-wave frequencies that is used to 
infer snowfall rates over land. The retrieved snowfall rates are qualitatively validated 
against ground-based radar observations of the same March 5-6, 2001 New England 
Blizzard. 
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Abstract: 

Because microwave brightness temperatures emitted by snow covered surfaces are highly variable, 
snowfall above such surfaces is difficult to observe using window channels that occur at low 
frequencies (v < 100 GHz). Furthermore, at frequencies v < 37 GHz, sensitivity to liquid 
hydrometeors is dominant. These problems are mitigated at high frequencies (v > 100 GHz) where 
water vapor screens the surface emission and sensitivity to frozen hydrometeors is significant. 
However the scattering effect of snowfall in the atmosphere at those higher frequencies is also 
impacted by water vapor in the upper atmosphere. This work describes the methodology and results of 
physically-based retrievals of snow falling over land surfaces. 

The theory of scattering by randomly oriented dry snow particles at high microwave frequencies 
appears to be better described by regarding snow as a concatenation of “equivalent” ice spheres rather 
than as a sphere with the effective dielectric constant of an air-ice mixture. An equivalent sphere snow 
scattering model was validated against high frequency attenuation measurements. Satellite-based high 
frequency observations from an Advanced Microwave Sounding Unit (AMSU-B) instrument during 
the March 5-6, 2001 New England blizzard were used to retrieve snowfall over land. Vertical 
distributions of snow, temperature and relative humidity profiles were derived from the Pennsylvania 
State University-National Center for Atmospheric Research (PSU-NCAR) fifth-generation Mesoscale 
Model (MM5). Those data were applied and modified in a radiative transfer model that derived 
brightness temperatures consistent with the AMSU-B observations. The retrieved snowfall distribution 
was validated with radar reflectivity measurements obtained from the National Oceanic and 
Atmospheric Administration (NOAA) National Weather Service (NWS) ground-based radar network. 
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1. Introduction: 


Measurement of global precipitation is one of the goals of climate studies. Although most global 
precipitation occurs as rainfall, snowfall plays a significant role in the extra-tropical hydrological 
cycle. Snow, falling early in winter, can retard freezing of the underlying soil thereby allowing 
subsequent melt water to penetrate the ground. Conversely, if the ground freezes because snow falls 
late in winter, flooding may ensue from run-off during the spring thaw. Snow also serves as a reservoir 
of water that can be released later in the year to support agriculture and hydroelectric power 
generation. Snowstorms can also be hazardous for transportation and other economic activities. One of 
the most important challenges for future satellites is to detect snowstorms from space. This study 
presents a physical model of radiation at millimeter-wave frequencies that seeks to infer snowfall rates 
over land by taking advantage of water vapor screening to obscure the underlying snow-covered 
surface. 

Mid-latitude precipitation occurs in a wide variety of forms from snow to drizzle to rain and to hail 
[1]. Storm types include multicell and supercell thunderstorms, extratropical cyclones and hurricanes. 
The well developed techniques to study tropical precipitation, using frequencies < 90 GHz, addressed 
rain occurring in nearly moist adiabatic environments. Extratropical cyclones are a completely 
different setting for precipitation than any type of tropical storm. Broad zones of frontal lifting 
produce broad sheets of clouds and precipitation that are mostly though not completely stratiform. 
These stratiform precipitation regions are similar microphysically to the stratiform regions of tropical 
mesoscale convective systems. However the generally cooler conditions often produce a melting layer 
that is near the earth’s surface. Under the coldest conditions, the precipitation reaches the surface as 
snow. 

Most spacebome remote sensing of snow has addressed the measurement of snow accumulation on the 
ground, see [2]. Snow within the atmosphere has only been derived from microwave radiometry over 
oceanic regions where the measurements were not affected by snow accumulated on the ground [3 J- 
[6]. Furthermore, most of the snow considered in the above studies are frozen particles above the 
melting layer or anvil ice clouds, not snow falling at the surface. Because snow accumulation on land 
affects the emission properties of the surface, the measurement of snowfall within the atmosphere has 
been difficult with radiometers that operate at frequencies less than 100 GHz where the atmosphere is 
relatively transparent and the sensitivity to frozen particles is lower than at higher frequencies. Snow 
falling over land has been derived from the brightness temperatures at frequencies where absorption 
occurs using empirical relationships by Kongoli et al. [7] and by Chen and Staelin [8]. Although such 
empirical relationships are operationally useful, physical models are needed to understand how the 
retrieved snowfall depends on the various ground and atmospheric factors that affect the measured 
brightness temperatures. To our knowledge, this is the first retrieval of snow falling over land based on 
a physical model. 

The Advanced Microwave Sounding Units (AMSU-B) radiometers on the NOAA 15, 16, 17 
spacecraft [9] have the channel set and resolution to resolve locally intense precipitation. The AMSU- 
B have 10 x ~25 km footprints and provide observations at 89, 150 and 183 ±1, ±3, +7 GHz. These 
channels are sensitive to both the water vapor (for surface screening) and the snow particles. The 
AMSU-B radiometer on NOAA- 15 initially encountered radio frequency interference from Qn-board 
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transmitters that were ultimately shut down in the autumn of 1999 so that reliable spacebome data at 
frequencies greater 100 GHz were available by the winter of 2000-01. This study presents a physical 
model that was used to derive snowfall over land from AMSU-B observations. 

2. Case Study: 

The blizzard of 5-6 March 2001 presented a unique opportunity to observe intense snowfall over land. 
That blizzard was one of the more intense snow storms of the season, depositing on the order of 50 cm 
of snow on much of Vermont, New Hampshire and northeastern New York with several stations 
reporting that 75 cm were deposited for the day. Both the NOAA 15 and 16 satellites observed this 
blizzard (NOAA- 17 was launched in June 2002). However, the best spatial and temporal coverage 
between available ground radar data and AMSU-B data was at 2300 UTC with the NOAA- 15 AMSU- 
B observations. 

2. 1 Radar Data: 

Fig. la shows a composite of the National Weather Service (NWS) operational weather radar 
reflectivity, Zeff(mm 6 /m 3 ) obtained from several ground stations over the Northeastern United States 
on 5 March 2001 at 23:00 UTC. Note that the limited range of the radar data does not extend far over 
the ocean area (Fig. la). The snowfall was greatest over Connecticut, Maine, Vermont, and New 
Hampshire. This composite of Z e ffis based on whichever of the lowest four antenna elevations yield 
the highest reflectivity. At ranges beyond 50 km from the radar, those elevation angles are usually 
0.5°. The heights at which those reflectivities are measured varies with distance from the particular 
radar, falling between 0.5 and ~ 2.5 km. Although the NWS operational radar data have well known 
limitations, in the absence of a preplanned field observation campaign, they provide readily available 
observations to compare to snowfall derived from microwave brightness temperatures. 

The radar reflectivity data were smoothed with a 16 x 16 km template to match the finest spatial 
resolution of the AMSU-B channels. The center points of the smoothed radar data matched those of 
the AMSU-B latitude and longitude center points for each footprint. Since the NWS radar reflectivity 
resolution is very fine however its latitude and longitude mapping was not precise (offsets by no more 
than 0.1 degrees). Averaging the NWS image to the AMSU-B resolution tended to smooth any effects 
of location mismatch. The maximum reflectivity in the smoothed radar reflectivity data over the land 
is ~37 dBZ. Depending on the relationships used, this reflectivity can correspond to rain fall rates of 
up to 6.7 mm hr* 1 (or snowfall rates of 67 mm hr' 1 ). 

2.2 Microwave Data from AMSU- B: 

Brightness temperatures, Tbiso, measured by the 150 GHz channel of the NOAA-15 AMSU-B, at 2302 
UTC 5 March 2001, are shown in Fig. lb. Note the cold brightness temperatures (< 240K) in the 
blizzard region over Vermont and New Hampshire (near 42-44N, 71-74W), and the absence of 
contrast in surface features such as the Gulf of St. Lawrence (near 47N, 72W) and the Great Lakes 
(near 44N, 77W). This brightness temperature distribution is similar to the radar Z e ff values shown in 
Fig la. Figs, lc, d, e, and f show the distributions of the 89 and 183 + 7, +3, +1 GHz brightness 
temperatures, (Tb 89 , Tbi 83 + 7 , Tb| 8 3 ± 3 and Tbi 8 3 ±i, respectively) measured from AMSU-B at 2302 UTC 
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on 5 March 2001. It can be seen in Fig. lc that some of the cold Tb 8 9 values scattered over Canada may 
have been caused by accumulated antecedent snow. At 89 GHz, it is difficult to distinguish snow in 
the atmosphere from snow on the ground. Note that the Great Lakes and the St. Lawrence River, that 
are evident in Fig. lc, are screened by water vapor in the 183 GHz channels in Figs. Id, e, f. Although 
brightness temperatures are more affected by the water vapor profile, evidence of the snow still 
persists in the 183 + 3 GHz channel. The Tbi 83 ±i measurement reveals little indication of the blizzard, 
being mainly responsive to water vapor at altitudes above the heaviest precipitation. It should be 
noted that the core of the snowstorm is at about 35° from nadir in the AMSU-B scanning pattern. 

3. Snow Retrieval Method: 

This study seeks to derive characteristics of snow whose electromagnetic properties are consistent with 
microwave brightness temperatures at several frequencies provided by the AMSU-B sensors. 
Brightness temperatures are computed from an Eddington approximation of the second kind. That 
radiative transfer model employs information generated from a six parameter model of the atmosphere 
associated with snow storms. Three of the six parameters are allowed to vary to generate different 
snow cloud and surface conditions, the other three parameters are set to fixed values based on statistics 
from a cloud resolving model and external measurements. It is expected that there are only 3-4 degrees 
of freedom in the five AMSU-B brightness temperature channels. 

The first and second of the three variable parameters are used to adjust the distributions of the vertical 
structure of snow mass and relative humidity and are generated from the Pennsylvania State 
University-National Center for Atmospheric Research (PSU-NCAR) fifth-generation Mesoscale 
Model (MM5). During the retrieval the parameterized snow model profile mass is allowed to extend 
beyond the confines of the MM5 simulations. This is necessary because the maximum surface snow 
mass over land produced by MM5 was only 0.7 g m‘ 3 corresponding to a lower reflectivity of 30 dBZ 
than seen in Fig. la. Thus the retrieved profiles are not just weighted combinations of the MM5 
profiles and they may vary more widely than the MM5 model profiles. The microphysical model of 
the falling snow also rests upon independent millimeter- wave attenuation measurements. The third 
adjustable parameter defines the fraction of snow ground cover and generates a composite surface 
emissivity using previously reported emissivities of snow and soil covered surfaces. 

The three fixed parameters include height levels, temperature profiles, and snow size. These 
parameters, along with the variable parameters, produce snow cloud characteristics used to generate a 
database of brightness temperatures that would be observed at the AMSU-B frequencies using forward 
radiative transfer calculations. The optimal estimate of the snow parameters is derived from the best 
match between computed and measured brightness temperatures at all AMSU-B frequencies. 

Snowfall mass distributions are derived by matching measured brightness temperatures to those 
derived from radiative transfer models applied to MM5 temperature and adjusted humidity profiles. 
The snowfall in the retrieval model is represented by equivalent spheres whose diameters are of the 
same order of magnitude as those needed to model attenuation measured during previous field 
experiments. Finally the retrieved snowfall is compared to the radar reflectivity measurements 
provided by the NWS operational radars to estimate the validity of the retrieval. 
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3. 1 MM5 Mesoscale Model: 


The MM5 model described by [10], [1 1] was used in this study. A description of the MM5 model used 
here was presented by Chang et al. [12]. The model domain was centered at 35°N x 70°W; it consisted 
of 100 x 120 grid points with 40 km separation. The simulation generated profiles of temperature, 
relative humidity, cloud water, rain, cloud ice, snow, hail, and graupel mass at each of the grid points 
in the volume of the model domain. 

The MM5 model was initialized at 0000 UTC 5 March 2001, and model integration was performed for 
a period of 48 h. The initial and boundary conditions for the experiment were obtained by interpolating 
the National Center for Environmental Prediction (NCEP) 2.5° x 2.5° global analysis for 0000 UTC 5 
March 2001 to the model grid points and then enhancing those with coincident rawinsonde soundings. 
The sea surface temperature data were also obtained from NCEP global analysis and were kept 
unchanged during the model integration. However, the land surface temperature was predicted using 
surface energy budget equations in which the effect of short and long wave radiation and cloud 
radiation were included. 

Any rain erroneously generated over coastal regions by the MM5 simulation was assumed to be snow 
because all of the reported temperatures in those areas were below - 2° C. The National Lightning 
Detection Network (NLDN) [13] observed some lightning off-shore, but not over New England. The 
cold land surface temperatures appeared unable to produce sufficient convection to produce lightning 
and perhaps graupel. Although not shown, the shape (but not the location) of the spatial distribution 
over land and ocean from the MM5 model compares well with that of the radar reflectivity, Z e fr, 
shown in Fig. la. 

3.2 Snow Cloud Parameterization: 

Three fixed parameters in the retrieval include height levels, temperature profile, and snow size. The 
height levels and temperature profiles were set to fixed values provided in Table I. The temperature 
profiles were determined by the average of the MM5 temperature profiles over land. The variability of 
the MM5 temperature profiles over land is minimal (about 10 K at lower altitudes). The three variable 
parameters and the fixed snow size will be described in detail in the following sections. 

3.2.1 Relative Humidity. Parameters 

Profiles of the relative humidity with respect to ice used in the retrieval were modified by adding a 
scaled addition to a minimum relative humidity profile generated by the MM5 over land. Those 
profiles were interpolated between the values cited in Table 1 to yield: 

RH(z) = RH m i n (z) + ARH(z) r (1) 

where ARH(z) = RH max (z) - RH m j n (z) is the maximum range of acceptable relative humidity values 
from MM5 and z denotes the height. The adjustable relative humidity scaling parameter r ranges 
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from 0.0 to 1 .0 in 0. 1 unit increments. The value of r was determined from the retrieval by minimizing 
the difference between observed and retrieval calculated brightness temperatures 

Table I. Atmospheric parameters inferred from the MM5 model of the March 5-6, 2001 blizzard. 


Height 

(km) 

T 

(K) 

RH (%) 
RH mi „ ARH 

M(z), 

Normalized 
Snow Profile 

0.02 

267.50 

80 

20 

1.00 

0.5 

267.13 

70 

30 

0.95 

1.0 

266.75 

60 

40 

0.90 

2.0 

266.27 

20 

80 

0.76 

3.0 

265.23 

11 

89 

0.61 

4.0 

261.72 

9 

91 

0.51 

5.0 

255.77 

6 

94 

0.35 

6.0 

^248.64 

4 

96 

0.20 

8.0 

232.51 

2 

98 

0.06 

10.0 

221.37 

2 

98 

0 

12.0 

222.52 

2 

34 

0 

14.0 

220.94 

2 

18 

0 

16.0 

216.61 

2 

16 

0 


3.2.2 Surface Emissivitv Parameter 

The radiative transfer equation also requires knowledge of the emissivity of the variable surface 
features including accumulated snowfall. The boundary conditions were determined partially by the 
accumulated antecedent snow whose emissivities, £ s , for deep dry snow at the relevant frequencies are 
obtained from [14]. Although [14] only measured emissivities to 150 GHz, they did provide curve fits 
for deep dry snow extending to 200 GHz. The 1 83 GHz emissivity was extracted from these curve fits. 
The values of e* at a 53° viewing angle are 0.64, 0.724, 0.8 for 89, 150, and 183 GHz respectively [14]. 

The emissivity used in the radiative transfer model is a weighted mean of the emissivity of snow cover, 
Es, and that of £<,. The £o is an average of emissivities from bare soil, frozen soil, and winter 
forest/conifer and is 0.98 for all frequencies [15]. The effective emissivity is thus 

£ = ffEs + O-ffjEo (2) 

where ff is the fraction of the ground covered by snow and has 6 values: 0.0, 0.2, 0.4, 0.6, 0.8, and 1.0. 

The ff parameter mostly affects the 89 GHz brightness temperatures that permit the earth’s surface to 
be seen through light precipitation. 
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3.2.3 Snow Mass Parameter 


Representative snow mass profiles were taken from the MM5 snowing profiles over land and 
normalized with respect to their surface snow mass, to obtain a normalized profile, M(z), shown in 
Table 1. Note that we consider the surface snow mass to be the snow mass at 20 m above the surface 

(see Table 1). Since the MM5 snow profiles truncated to 0 g m' 3 at a 10 km height, the M(z) profile 
has 0 mass at a 10 km height (See Table 1). Snow mass profiles used in the subsequent retrieval are 
scaled by a factor m in g m' 3 . The snow mass scaling parameter m could assume 39 values: 0.0, 0.02, 
0.065, 0.1, 0.2, 0.4, 0.6, ..., 6.8, 7.0 g m' 3 . These values convert to surface melted snowfall rates 
ranging from 0 to ~32 mm hr' 1 , or equivalently, to reflectivities ranging from 0 to ~48 dBZ. These 
snowfall rates, although beyond the surface maximum of 0.8 g m' 3 given by the MM5 model, provide 
enough variability for the range of values seen in Fig. la. The snow mass content used in the retrievals 
is 


M s (z) = m M(z) . 


(3) 


3.2.4 Snow Size Panimeter Selection 

The greatest challenge of snowfall studies is determining the electromagnetic properties of the wide 
variety of shapes and sizes of snowflakes. Microwave scattering and emission properties of 
hydrometeors depend on their size, shape, density, and whether they are frozen, melting or liquid. Wet 
snow and sleet were reported along the New England coast, but the mean temperatures encountered in 
New Hampshire and Vermont remained around -5 C, and reported maxima were only -2 C, we 
therefore disregarded melting over inland regions throughout the entire day. Furthermore, our case 
study observations are at 2300 UTC or 1900 local time when the temperatures are less than the 
daytime highs. Snow particles in this study were therefore regarded as dry containing no melt water. 

We first assumed that the mass densities of the MM5 model consisted of large fluffy snow particle size 
distributions (i.e., 10% ice, 90% air) as suggested by Rutledge and Hobbs [16]. Effective medium 
mixing theories have been used by [5], [6], [17], [18] to represent the dielectric constant of snow at 
frequencies below 90 GHz [19]. Using these models for the higher frequencies seemed to produce 
inappropriate electromagnetic characteristics. For example, the low density particles represented by 
effective medium mixing models did not provide enough scattering because the asymmetry factor was 
too large. Large asymmetry factors increased forward scattering of the radiation from the warm lower 
layers so that computed brightness temperatures were too warm. Although the finite difference time 
domain method [20] can be used to compute the scattering characteristics of non-spherical particles at 
any frequency the shape of the frozen crystal habit can only be crudely estimated, so that a simpler 
approach appears to be justified. 

One such simpler approach is the procedure of Grenfell and Warren [21]. Grenfell and Warren (G-W) 
represented randomly oriented frozen particles as a concatenation of “equivalent” ice spheres whose 
effective diameter was determined by the ratio of the volume-to-surface-area, (V/A). 


D eff = 6V/A 
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( 4 ) 



The G-W procedure transforms inhomogeneous non-spherical (e.g., fluffy) ice particles into an 
ensemble of solid ice spheres; this greatly simplifies the determination of the scattering properties of 
the hydrometeors. Representing irregular particles as ensembles of equivalent spheres, has been used 
in the cirrus cloud infrared radiation modeling community and their properties have been analyzed by 
[22]-[23]. G-W demonstrated that equivalent spheres can adequately describe the transmittance and 
reflectance of diffuse infrared radiation through randomly oriented prisms. Moreover, [24] showed 
that the equivalent sphere approach accounts for the shape transition from needles to plates. It is 
noteworthy that such effective diameters are mainly determined by the small dimensions i.e., the 
thickness of large disks or the diameters of long cylinders (Eq. 3. in G-W) rather than the maximum 
dimensions that are most frequently measured. 

The distributions of effective diameters, D e jf (mm) based on the G-W model were represented by a 
Gamma function of order 1. 


N(D e ^ - No D e jf ex/?(-AcwA#). 


(5) 


The No is related to the G-W slope parameter Aow and the snow mass density, M s (g m' 3 ) at each 
height, z: 


= (D„)D a 'dD a 

0 0 

where pow = pi = 0.917 g m' 3 so that 

XT _ M s {z)\q W 
■‘ V 0 ” A 

^Pgw 


( 6 ) 

(7) 


Note that the No value varies with differing snow masses at each height, through Eqn. 7 and that the 
number of “equivalent” spheres is greater than the actual number of non-spherical scattering particles. 
The effective polydispersion diameter weighted over such a size distribution is: 

tf>„> (8) 

J 'GW 


Once the G-W size distribution was specified, the effective diameter, < D eff > was determined by 
employing measurements obtained from other snow events to infer the attenuation coefficient (and 
hence < D e ff>) of snow to use in the electromagnetic scattering model. The size parameter < A#> was 
derived by applying Mie theory to compute the extinction coefficient, ksext (km 1 ) of snow as a function 
of frequency, v(GHz). The attenuation per mass is: 

A( v, M s ) = 10 log l0 (e) ksext / M s , (9) 
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Results of those calculations were compared to A(V, M s ) measured for several winter seasons for 
frequencies between 96 and 225 GHz by [25], [26]. Fig. 2 compares the attenuation per mass as a 
function of frequency measured during the samples with curves generated by Mie theory for <D e ff > 
values of 0. 1mm and 0.06mm. It is well known that the crystal habit varies with height and that <D e jj > 
diminishes with height [4]. A simplification was invoked because it is difficult to justify the 
introduction of diverse shapes in the model in the absence of measurements. Thus snow particles were 
assumed to be spherical and to be represented by: 


<D e ff> = 0. 10 mm for 0.0 < z < 0.5 km 

= 0.06 mm for 0.5 <z< 10.0 km. (10) 

These dimensions are characteristic of the small crystal dimensions, i.e. the thickness of plates or 
diameters of needles [27], [28]. Furthermore, the <D e ff> = 0.1 mm below 0.5 km better matched the 
ground-based observations of Nemarich et al. [25] and Wallace [26], while the <A# > - 0.06 mm 
above 0.5 km yielded brightness temperatures that best match the AMSU-B observations. In fact, 
section 4.2 will show that the 150 and 183 + X GHz brightness temperatures are especially sensitive to 
the particles found at higher altitudes. 

3.3 Radiative Transfer Model: 

The radiative transfer calculations are an integral part of the retrieval method since they are used to 
minimize the error between the observed and the calculated brightness temperatures. The radiative 
transfer model requires vertical profiles of pressure, temperature and humidity as well as cloud water 
and cloud ice and precipitation. These quantities were obtained from the snow cloud parameterization 
and they were introduced into the Mie theory model. Although the scattering by snow is nearly 
conservative, the actual albedos for single scattering that appear in the radiative transfer model were 
not, because those were determined by absorption mainly due to water vapor. The albedos for single 
scattering thus varied with height as the constituents varied. 

The radiative transfer model used to compute brightness temperatures given a hydrometeor profile is 
the delta-Eddington model. Microwave radiances 35° off nadir were computed from a second order 
Eddington approximation [29] with delta scaling [30] for plane-parallel clouds. This permitted the 
radiances to be calculated analytically by representing them by a series of Legendre polynomials 
truncated after the first order. Those radiances were then inserted into the source function of the 
transfer equation to compute the brightness temperatures. Smith et al. [31] showed that such an 
algorithm produced acceptable brightness temperatures at a 53.1° viewing angle, but that it yielded 
unacceptable errors in brightness temperatures viewed at nadir. Accordingly the phase function was 
approximated by representing the forward scattering component with a Dirac delta function in addition 
to a constant term and a term proportional to the cosine of the scattering angle. That procedure enabled 
both first and second moments of the approximate phase function to match those of a Henyey- 
Greenstein phase function. This model is identical to the previously cited Eddington second-order 
approximation except that the profiles of asymmetry factor, extinction coefficient and albedo for 
single scattering were scaled by transformed parameters. Kim et al. [32] and [31] showed that this 
transformation reduced the computational errors for both nadir and 53.1° viewing angles. We therefore 
assume that the radiative transfer model is valid for AMSU-B scanning angles between those limits. 
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3.4 Snowfall Retrieval: 

The parameters r, ff and m were found by an optimization that sought the minimum of: 

T^r, f, m) = £ [Tbj- Tbj 0 ] 2 = minimum (11) 

where Tbj and Tbj° are the computed and observed brightness temperatures respectively, and the 
summation is over the five AMSU-B frequency channels. Since, for the current retrieval algorithm 
configuration, there are 11 r, 6 ff, and 39 m values, the minimization occurs over the computed 
brightness temperature vectors associated with these 2574 potential cloud profiles. When additional 
information is acquired about the snow or other aspects of the storm conditions, these can be included 
in Eq. (1 1) to further constrain the optimization. Furthermore, future versions of the retrieval algorithm 
will be modified for iterative optimization to allow any values of r, f, m to be selected solely by Eq. 
( 11 ). 

4. Retrieval Results and Analysis: 

The retrieval was performed for the March 5-6, 2001 blizzard. The land surface temperature was 
assumed to be 267.5 K throughout the entire scene and the surface pressure was assumed to be 
lOlOmb. This surface temperature is a few degrees cooler than the reported inland daytime highs of - 
2° C, since the local time of the AMSU-B observations was shortly after sunset. Figs. 3 a, b, c, and d 
illustrate the spatial distribution of brightness temperature errors that contributed to the ^(r, ff, m) 
residuals at 89, 150 and 183 + 7, 3 GHz respectively. Because this study is confined to the 
determination of snowfall over land, results over water are grayed out. The theoretical and 
experimental brightness temperatures agreed within approximately + 5 K over most pixels at all of the 
AMSU-B frequencies; including 183 + 1 GHz which was not shown. The significance of the retrieved 
variables will be discussed in the following sections. 

4. 1 Retrieved Parameters: 

Fig. 4a shows the retrieved snowfall mass density near the surface, M i (0.02km) = m M(0.02km). It is 
evident that the spatial distribution of the snowfall mass is qualitatively similar to the radar 
reflectivity displayed in Fig. la. Figure 4b shows the distribution of the relative humidity at a 1 km 
altitude. The retrieved snow mass and relative humidity distributions are somewhat noisy. If more 
information were known about the profile (such as through vertical radar profiles and in-situ 
measurements) it is expected that these variations would be reduced. Potential sources of this noise 
are discussed in Section 4.3. 

Figure 4c shows the distribution of the parameter f, the fraction of snow cover on the ground within 
the AMSU-B field of view. Water in the Montreal, Quebec, and the St. Lawrence River regions (near 
45-46N, 73-76W) may skew the retrieved snow cover fractions because the water surfaces have a 
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lower emissivity than dry ground surfaces assumed in Eq. (2). The retrieval algorithm erroneously 
compensates for these lower emissivities (that produce cooler Tb) by increasing the snow fraction. 

4.2 Weighting Functions: 

Temperature weighting vector profiles [33] can be used to understand contributions to brightness 
temperatures from each layer of the atmosphere, the cosmic background, and the ground. The 
brightness temperature at each frequency is the integrated sum over all heights of the product of the 
weighting vector value and the atmospheric temperature plus the contributions from the ground and 
cosmic background temperatures: 


Tb = ]T(z)W(z)dz = T 0 W 0 + T C ,W CI + £ TWw. ( t)tu , 

0 *-l 


( 12 ) 


where the subscripts 0 (at z - 0), and CB denote ground surface and cosmic background, T(i), W A u(i) 
denote the temperature and weighting vector value for level / of the cloud profile that consists of I 
levels, and A z is the height increment between level i and level /-l. The weighting vectors depend on 
the profiles of the atmospheric cloud constituents and will change for various representative profiles. 
The weighting vectors Wall, Ws, an d Wcb are defined in [33] where Wall is denoted W (n) . The Wall 
includes the effects of multiple scattering by snow (S), and absorption by snow, water vapor (WV), 
nitrogen (N 2 ), and oxygen (O 2 ). 

In order to determine the contributions to the brightness temperature value from snow, water vapor, 
and other constituents, W A llT(z) (with units of K km' 1 ) at each height is revised: 


W ALL (z)T(z) = M s (z)W SNOW (z) + RH(z)W wv (z) + W olher (z)T(z ) . (13) 

We then define Ws, Wwv, and W 0 tt, er as: 


where 


jy ( z ) — ^Sext (z)TXz)F^(z) 

k ex,( Z ) M s(z) 


W^(Z) 


^(z)T(z)IT m (z) 

k ex ,(z)RH(z) 


Jir / \ _ [^02 ( z ) + ktf 2 (^}WaLL ( Z ) 

n other \ 2 )- . f v 

*«(*) 


(Kkm' 1 m 3 g _1 ) 
(Kkm" 1 %'’) 
(km' 1 ) 


(14) 

(15) 
(16) 


Kxt - ks.ext + kwv + ko2 + k N 2 . (17) 

The ks,ext term is the sum of the absorption and scattering coefficients for snow found in Eq. 9, kwv is 
the absorption coefficient for water vapor, and ko 2 and km are the absorption coefficients for O 2 , and 
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N 2 . Note that k ex , is also used in integrals within Wall as the atmospheric opacity. This means that Ws 
is not entirely free of the effects of the water vapor and other constituents, and likewise for W W v and 
W 0 ,her however, they do respond to most of the information about the snow, water vapor, and other 
constituents, respectively. 

In order to illustrate the physics of this retrieval, two representative retrieved profiles are studied. The 
first retrieved profile is extracted near the heaviest snowfall region, at 42.52N, 72.036W. The second 
profile at 40.77N, 72.36W is a lighter snowfall case and it is also much drier. For the first profile, Eq. 
(1 1) selects r = 0.7, m = 2.6 g m* 3 , and f = 0.8 and for the second profile r = 0.3, m = 0.6 g m'\ and f 
= 0.4. The retrieved snow and relative humidity profiles and surface emissivity can be computed using 
Eqns. (1), (2), (3) and the profile information in Table 1. 

Table II compares 4he observed and computed brightness temperatures for these two representative 
retrieved profiles. This table shows fairly good agreement for all but the profile 2 183 ± 1 GHz 
channel where there is a 5 K difference. This disagreement is likely due to a mismatch in the relative 
humidity for the second profile. 

Table II. Brightness temperature values (in K) for two representative profiles. 



Profile 1 

Profile 1 




AMSU-B 

Retrieved 



89 GHz 

209.2 

206.5 

233.9 

232.0 

150 GHz 

185.5 

185.7 

221.4 

219.5 

183 ± 1 GHz 

236.8 

237.2 

241.4 

246.3 

183 ±3 GHz 

234.1 

232.9 

244.3 

247.4 

183 ±7 GHz 

210.1 

209.4 

235.1 

236.0 


Plots of Wall(z)T(z ), Ws(z)M s (z), and Wwv(z)RH(z) are shown in Figure 5 for the two representative 
profiles. These plots can be used to determine the brightness temperature (by integrating over the 
heights and adding the cosmic background and surface contributions provided in Table 3). The 
Wwv(z)RH(z) plots (Fig. 5c and f) show that the contribution of water vapor to the 183.3 ±1 and ±3 
GHz brightness temperature is high at higher altitudes (above 5km). For the other frequencies, the 
water vapor contribution is inversely correlated with the amount of snow in the cloud. By examining 
Fig 5b and d, it is noted that a decrease in m from 2.6 (Fig. 5b) to 0.6 g m' 3 (Fig. 5d) reduces the 
contribution from snow by about one half. Furthermore it should be noted that the 1 83 GHz channels 
response to snow above the Earth’s surface, while the 89 and 150 GHz respond more intensely near 
the surface. The 89 and 150 GHz channels are also more sensitive to the surface emissivity, so some 
tradeoff must be made between the surface contributions and the snow contributions and 
contamination of the signal may occur. The 150 GHz channel is a little more immune to surface 
effects if the snow mass profile is large enough. 
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Table III. Surface and cosmic background weights and brightness temperature contributions for three 

retrieved profiles. 



Profile 1 

Profile 1 

Profile 2 

Profile 2 


Bn 

TcbWcb 

ToWo 

TcbWcb 



0.5 

180.8 

0.37 

150 GHz 

25.7 

0.8 

110.9 

0.5 

183 ± 1 GHz 

~ 0.0 


0.01 

0.06 

183 ±3 GHz 

0.02 

0.29 

0.65 

0.15 

183 ±7 GHz 

1.4 

0.55 

16.4 

0.33 


The ToWo in Table III show that the surface contribution from the 89 GHz channel is larger than that of 
the other channels and increases as we move from profile 1 to profile 2. In fact, the surface contributes 
at least half of the value of the resulting brightness temperatures for 89 GHz. A similar effect is seen 
for the 150 GHz surface weighting values, though heavy snow will obscure the surface for 150 GHz 
(Profile 1). The 183 ± X GHz channels receive very little of their resultant brightness temperature 
values from the surface. All of the cosmic background contributions, TcbWcb, are small. 

4.3 Radar Reflectivity vs. Melted Snowfall Rain Rate 

A pixel area matching technique similar to that described by [34] was used to relate the radar 
reflectivity, Zeff (mnr m' 3 ), derived from the dBZ e fr over land shown in Fig. la to the retrieved 
snowfall mass distribution, A/*(0.02km) retrieved from the microwave data shown in Fig. 4a. The 
procedure selected a number of pixels that exceeded a given snowfall mass and the same number of 
pixels that exceeded a particular radar reflectivity. Threshold values of each of these quantities were 
tabulated. The terminal velocity of snowflakes was assumed to be ~ 1 ms' 1 so that the melted snowfall 
rates, R (mm h' 1 ), could be derived by multiplying A/ V (0.02km) by that terminal velocity. 

Fig. 6 presents dB Z e ff as a function of log(/?) derived from this pixel matching technique for the 
retrieval results reported in Fig. 4. These results are compared to some representative Z e fr -R 
relationships [35]-[38] that showed that the coefficients of snow Zefr -R relationships could vary over a 
wide range. The coefficients in such Z e ff -R relationships have been refined over the years as snow 
measuring techniques have become more sophisticated. The comparison of the retrieved relationship 
to some previously published relationships is surprisingly good. While this is not a direct validation of 
the retrieval results, it does show that this physical model enables retrievals to fall within existing 
measured and empirical relationship bounds. 

4.4 Error Discussion: 

Because there are only 3-4 degrees of freedom associated with the AMSU-B channels and few 
additional measurements to further constrain the retrievals, the number of parameters retrieved is 
limited. With a limited retrieved parameter set, assumptions were made about the vertical structure of 
the snow cloud, including fixed vertical profiles of temperature, water vapor, and snow mass shape. 
Unfortunately, these assumed profiles, taken from the MM5 cloud model simulations, could not be 
validated by in situ observations. Further constraints were made on the allowed variability of the snow 
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particle characteristics, by keeping the snow particle effective diameter relatively fixed throughout the 
profiles (but allowing the snow mass index, m to vary). 

Determining the electromagnetic (scattering, absorption, and asymmetry) properties of the snow leads 
to additional assumptions. The ones used herein include assuming dry snow, that the Grenfell and 
Warren [21] approach to redistributing low-density fluffy snow into small ice spheres is appropriate 
and that the snow attenuation measured by [25] and [26] applied to the present case with a minor 
modification of the particle size. Fig. 2 showed that the attenuation of the particles assuming the G-W 
approach is relatively consistent with independently measured attenuation of snowfall. Other research 
[3, 39] on anvil ice particles and snowfall has shown that the retrieved particles sizes can be smaller 
than expected giving credence to the applicability of the Grenfell and Warren approach. A detailed 
future analysis should include in situ sampling of snow particles that simultaneously determines size 
distributions and electromagnetic characteristics. 

The noise in the retrieved snow mass, relative humidity indices and surface emissivity may be 
attributed to some of the assumptions made in this physically-based retrieval: 

» 

1 . Selection of a single shape, M(z), for the snow mass profiles. Eq (3). 

2. Selecting a step function to describe the snow size distribution with height, Eq. (10). 

3. Assuming that equivalent spheres suffice to describe the electromagnetic properties of these 
frozen hydrometeors. Eq. (4-7). 

4. Assuming that the snow measured by [25], [26] had the same habit as that at the surface in this 
New England blizzard. 

5. Ignoring frozen hydrometeors other than snow and assuming that snow was dry. 

6. Allowing the relative humidity index to vary independently of the snow mass index. 

7. Assuming that topography had no influence on the surface emissivity other than that represented 
by Eq. (2). 

Another source of error and uncertainty is the conversion of the NOAA NWS composite data into 
reflectivities associated with snow falling near the ground. Radar validation of the melted snowfall 
rates retrieved from microwave brightness temperatures are fraught with large uncertainties both in 
terms of calibration of the operational NWS radars and the choice of the appropriate height from 
which measurements were blended to produce the mosaic of reflectivity values. While, the 
reflectivities presented in summaries of operational NWS radars were not calibrated for research; they 
provide information about the intensity and horizontal structure of the storm. In the absence of a 
dedicated measurement campaign the snowfall melted rain rates yielded a Z e ff-R relationship that was 
consistent with previous measurements. 

5. Discussion and Conclusions: 

A physically based retrieval algorithm was developed to estimate snowfall over land. The retrieval 
algorithm relied on a multi-parameter cloud model to generate the vertical structure of a snow cloud, 
including snow mass, snow particle effective diameter, and water vapor. The MM5 cloud simulation 
was used to provide useful statistics for generating those cloud characteristics. Ground-based 
attenuation measurements were used to characterize the equivalent sphere snow particle size used 
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herein. The snow cloud profile and surface emissivity were then used in radiative transfer calculations 
that were optimized against AMSU-B observations at 89, 150 and 183. 3±7, ±3, and ±1 GHz. For each 
pixel in the image, the multi-parameter cloud parameterization that produced brightness temperatures 
that best fit the AMSU-B observations was selected as the retrieved profile. The present study 
demonstrated: 

1 . Microwave radiometric channels operating at frequencies greater than 89 GHz provide information 
on snowfall over variable land surfaces because the surface emissivity is screened by water vapor 
absorption at those frequencies. 

2. An electromagnetic scattering model of randomly oriented snow particles was adequately 
represented as equivalent spheres whose diameters were mainly determined by the small 
dimensions of the snow particles as suggested by the work of Grenfell and Warren [21]. That 
model accounted for measured values of attenuation per unit mass between 96 and 225 GHz. 
Inserting the G-W “equivalent” ice spheres in a delta-Eddington radiative transfer model yielded 
brightness temperatures at 89, 150 and 183 + 1, +3, +7 GHz that were consistent with values 
measured by AMSU within + ~5K. This was in part due to the fact that the diameters of the 
equivalent particles were small so that the asymmetry factor was also small (as might be expected 
from Rayleigh-like scattering). Small asymmetry factors reduce the transmission of snow layers, 
thereby achieving lower brightness temperatures than those produced by a low-density fluffy snow 
particle ice-air effective medium with larger asymmetry factors. 

3. Weighting vectors illustrated the relationships between the physical properties of the clouds (snow 
and water vapor characteristics) and the resulting brightness temperatures. 

4. Three variables used to adjust the snow mass, relative humidity and surface emissivity were 
sufficient to estimate snowfall rates consistent with NWS radar reflectivity measurements and to 
yield a Zefr -R relationship that was consistent with others reported in the literature. 

The number of retrieved parameters was kept to a minimum because there are only 3-4 degrees of 
freedom in the five AMSU-B channels. This study emphasizes the need for a dedicated set of 
coincident observations that include microwave as well as microphysics measurements. Field 
campaigns are needed to measure the high frequency electromagnetic properties of snow along with 
the habits of frozen hydrometeors to yield parameters that we were forced to derive from disparate 
observations. Such measurements need to include the small as well as the large dimensions of frozen 
hydrometeors. More realistic retrieval procedures can be developed when additional information 
becomes available. 
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Table Captions: 

Table I. Atmospheric parameters inferred from MM5 model of the 5-6 March 2001 blizzard. 

Table II. Brightness temperature values for three retrieved profiles. 

Table III. Surface and cosmic background weights and brightness temperature contributions for three 
retrieved profiles. 

Figure Captions: 

Fig. 1 (a) Radar reflectivity (dBZ) obtained from the NWS operational radar composite measured at 
23:00 UTC 5 March 2001 and AMSU-B brightness temperatures observed 23:02 UTC 5 March 2001 
(b) 150 GHz, (c) 89 GHz, (d) 183 + 7 GHz (e) 183 + 3 GHz, and (f) 183 + 1 GHz. 

Fig. 2 Comparison between the attenuation per unit mass density (dB/km/ gm/m 3 ) of snow versus 
frequency. Calculations (solid lines) are compared to various measurements by (o) Wallace, 1988 
(1981 measurements), (X) Wallace, 1988 (1982 measurements) and (□) Nemarich et al, 1988. 

Figs. 3 a), b), c), and d) Differences between optimized computed brightness temperatures and 
measured brightness temperatures at 89, 150, 183 + 7, and 183 + 3 GHz respectively. 

Figure 4 The distribution of retrieved (a) snow mass at 20 m above the surface, (b) relative humidity at 
1 km height, and (c) fractional snow cover. 

Figure 5: The total, snow, and relative humidity weighting vectors of two sampled retrieved profiles, 
(a, b, c) near the most intense snowfall region, and (d,e,f) near the storm boundary with low relative 
humidity. 

Figure 6 Measured NWS radar reflectivity, Z, versus melted snowfall rate, R, (mm h* 1 ) inferred from 

the M s (0) (~X~X~). Results are compared to Z-R relationships found by Sekhon and Srivastava ( 

), Vasiloff et al. ( - • - ) and by Fujiyoshi et al. ( ), and Boucher and Weiler ( ). 
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Fig. 2 Comparison between the attenuation per unit mass density (dB/km/ gm/m 3 ) of snow versus 
frequency. Calculations (solid lines) are compared to various measurements by (o) Wallace, 1988 
(1981 measurements), (X) Wallace, 1988 (1982 measurements) and (□) Nemarich et al, 1988. 
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Fig. 5 The total, snow, and relative humidity weighting vectors of two sampled retrieved profiles, (a, 
b, c) near the most intense snowfall region, and (d,e,f) near the storm boundary with low relative 
humidity. 
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Fig. 6 Measured NWS radar reflectivity, Z, versus melted snowfall rate, R, (mm h' 1 ) inferred from the 
M s (0) ( - X“X _ ). Results are compared to Z-R relationships found by [34]-[37]. 
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